Particle number fluctuations in a cloven trapped Bose gas at finite temperature 
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We study fluctuations in the atom number difference between two halves of a harmonically trapped 
Bose gas in three dimensions. We solve the problem analytically for non interacting atoms. In the 
interacting case we find an analytical solution in the Thomas-Fermi and high temperature limit in 
good agreement with classical field simulations. In the large system size limit, fluctuations in the 
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independently of the trap anisotropy. The occurrence of this maximum is due to an interference 
effect between the condensate and the non-condensed field. 
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I. INTRODUCTION 

Fluctuations appearing when counting the atoms in 
a given sub-volume of a quantum system, are a funda- 
mental feature determined by the interplay between the 
atomic interactions and quantum statistics. They can be 
used to investigate many-body properties of the system 
and in particular non-local properties of the g( 2 > pair- 
correlation function. These fluctuations were studied at 
zero temperature for quantumgases in different regimes 
and spatial dimensions in |l|-[5||. Sub-poissonian fluctu- 
ations appear for non-interacting fermions and for in- 
teracting bosons. A related issue in condensed matter 
physics is that of partition noise in electron systems @. 

In cold atoms experiments it is now possible to di- 
rectly measure the fluctuations in atom number within a 
given region, as done for example in Q for a quasi one- 
dimensional system. However finite temperature plays 
a major role in experiments. Very recently, experiments 
were done on an atom chip where a cold gas of Rb atoms, 
initially trapped in an single harmonic potential well, is 
split into two parts by raising a potential barrier. Accu- 
rate statistics of the particle number difference between 
the left and right wells Nl — Nr is then performed in 
the modified potential. By varying the initial tempera- 
ture of the sample across the transition for Bose-Einstein 
condensation, they observe a marked peak in the fluctu- 
ations of the particle number difference below the tran- 
sition temperature T c while shot noise fluctuations are 
recovered for T > T c . For T < T c they finally get sub- 
shot noise fluctuations due to the repulsive interactions 
between atoms Q. Here we show that the peak of fluc- 
tuations for T < T c is in fact a general feature already 
appearing in a single harmonic well if we look at the 
fluctuations of the atom number difference Nr, — Nr be- 
tween the left half and the right half of the trap along 
one axis. We show that, contrarily to what happens to 
the total number fluctuations, fluctuations in the parti- 
cle number difference Nl — Nr can be computed within 



the grand canonical ensemble without any pathology. In 
the first part of the paper we address the ideal gas case 
for which we find the complete analytical solution in the 
grand canonical ensemble. We derive the asymptotic be- 
haviors for T <C T c and T 3> T c and we explain the 
physical origin of the "bump" in fluctuations of the par- 
ticle number difference for T <T C . In the second part of 
our paper we then address the interacting case. 



II. IDEAL GAS: EXACT SOLUTION 

We consider an ideal gas of bosons in a three- 
dimensional harmonic potential. The signal we are in- 
terested in is the particle number difference Nl — Nr 
between the left and right halves of the harmonic poten- 
tial along one direction, as shown in FigfT] In terms of 
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FIG. 1: We consider fluctuations of the particle number dif- 
ference Nl — Nr between the left and the right halves of a 
three-dimensional harmonic potential. 

the atomic field operators: 

N L -N r = f - / 1>H ■ (1) 
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Due to the symmetry of the problem, Nl — Nr has a 
zero mean value. It is convenient to express its variance 
in terms of the unnormalized pair correlation function 



g^{v,v') = {^{v)^{v l )i>{v l )i:{v)) 



(2) 



Normally ordering the field operators with the help of 
the bosonic commutation relation pulls out a term equal 
to the mean total number of particles: 



Vav(N L -N R ) = (AO +2 
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(2), 



r,r') 



r£L Jr'eL 



I f 5 (2) (rV) 

JrEL Jr'eR 



(3) 



We assume that the system is in thermal equilibrium in 
the grand canonical ensemble with /3 = 1/ksT the in- 
verse temperature and \i the chemical potential. Since 
the density operator is Gaussian, we can use Wick's the- 
orem and express g^(r,r') in terms of the first-order 
coherence function g^'(r, r') = (pjft (r)-0(r')): 

( r , r ') = g M ( r , r')^ 1 ) (r', r) + gW (r, r)gV> (r', r') . 

(4) 

The#M(r,r') is a matrix element of the one-body density 
operator 



5 «(r,r') = (r'|— r 



1 



■ z -± e phi — i 1 ' 
where hi is the single particle Hamiltonian 

2m 
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To compute gi, it is convenient to expand the one-body 
density operator in powers of the fugacity z — e 13 ^ [9] : 



(r,r') = (r'|£ 



i.-iph 



z e 



Mr}- 



(7) 



On the other hand for a harmonic potential the matrix 
elements of e~P hl are known We then have: 



W(r,r') = E 



exp 



\2irh) 

mui a 
Ah 



3/2 



Yl [sinh(lr) a )] 



-1/2 



ot—x,y,z 



(,„ + r^tanhl% 



(8) 



where we introduced the geometric mean of the oscil- 
lation frequencies Co — (w I w J u z ) 1 '' 3 and r\ a — (3hu a . 
It is convenient to renormalize the fugacity introducing 
z = zcxp(— J2 a Vol/2) that spans the interval (0, 1). Af- 
ter some algebra [Tl|, the variance of Nl — Nr is ex- 
pressed as a double sum that we reorder as 



Var(Af L - N R ) = (N) 



c s z 



(9) 



with 



s_1 1 — ^ arctan J tanh^Zrfc) tanh \{s — l)rj x 



i=i 



n [i 



(10) 



with ci = 0. Correspondingly, the mean atom number is 
expressed as 



-i TT 1 + coth(;r? a /2) 



(11) 



1=1 a 



This constitutes our analytical solution of the problem in 
the grand canonical ensemble. 

In practice, the forms © and (fTTj) are difficult to han- 
dle in the degenerate regime, since the series converge 
very slowly when z — > 1. A useful exact rewriting is 
obtained by pulling out the asymptotic behaviors of the 
summands. For the signal we obtain the operational form 



Var(Af L - N R ) = (N) + Coo {N ) + ^(c s - Coo )z s (12) 
where 



= lim c s = 2 j 1 arctan \ / tanh ■ ^ 3 



i=i 



and (No) = z/(l — z) is the mean number of condensate 
particles. The mean atom number is rewritten as 



(N) = (N )+J2 



i=i 



-1 



coth(Z77 a /2) 



(14) 



In Fig. [2]we show an example of fluctuations of the parti- 
cle number difference for realistic parameters of an atom- 
chip experiment. It is apparent that fluctuations are 
weakly super-poissonian above T c and a marked peak of 
fluctuations occurs for T < T c . In the following sections 
we perform some approximations or transformations in 
order to obtain explicit formulas and get some physical 
insight. 



III. APPROXIMATE FORMULAS FOR k B T > hw 

In this section, we consider the limit of a large atom 
number and a high temperature fcgT 3> hw a for all a. 

Non-condensed regime: Taking the limit rj a <C 1 in Eq. 
(gll we get 



(N) 



k B T 

hcu 



9z{z) 



(15) 



where g a (z) = 'Yld=iZ l /l a is the Bose function. From 
this equation we recover the usual definition of the critical 
temperature T c : 



knT r 



[A7C(3)] 1/3 to, 



(16) 
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where £(3) = 33(1) with £ the Riemann function. Taking 
the same limit in © gives 



Var(JV L - N R ) ~ (N) 
At T = T c this leads to 



1 



g 2 (z) -g 3 (z) T 3 



Var(iV L - N R )(T C ) ~ A 



C(3) 
C(3) " 



1? 



1.37A. 



(17) 



(18) 



showing that the non-condensed gas is weakly super- 
poissonian, as already observed in Figj2] 

Alternatively, one may directly take the limit rj a — > 
in Eq.©, yielding 



1} (r,r')^E 



1=1 



;3/2 A 3 



n 



exp ■ 



Imul / r„ + r' N 



(19) 



This semiclassical approximation coincides with the 
widely used local density approximation, and allows to 
recover (fTB)l and (fTT]) . 

Bose- condensed regime: In this regime z — > 1 so we use 
the splittings (fT2"]) and (TT41 . Setting 1 = 1 and taking the 
limit rj a — > in each term of the sum over I in (|14[) , we 
obtain the usual condensate fraction 



(No) 
(N) 



T 3 

l -Tr 



(20) 



The same procedure may be applied to the sum over s in 
. The calculation of the small- 77 limit of Cqq requires a 




FIG. 2: (Color online) Normalized variance of the particle 
number difference Nl — Nr as a function of temperature in 
a cigar-shaped trap with uj v — ui z = 2uj x . The number of 
particles is (TV) = 6000 (black lines) and (TV) = 13000 (red 
lines). The inset is a magnification of the T > T c region. Solid 
lines: exact result ((9| together with (|11[) . Dashed line for T > 
T c : approximate result (|17[l together with (|15[l . Dashed line 
for T < T c : approximate result resulting from the improved 
estimates (|26|) and (|25|) . The temperature T is expressed in 
units of the critical temperature T c defined in Eq. (| lfcif) . 



different technique: Contrarily to the previous cases, the 
sum in (|13p is not dominated by values of the summation 
index <C 1 /r) a and explores high values of I ~ 1/t] x ■ As a 
remarkable conse quen ce, the local-density approximation 
fails in this case [l2j. We find that one rather has to 
replace the sum over I by an integral in (|13|) : 



+°° 21n2 
dlf(l) = 



(21) 



with f(x) = 2 — (8/71") arctan A k /tanh(x77 a; /2). This leads 
to the simple formula for T <T C : 



Var(7V L - N R ) ~ (A) 



-2 In 2 



C(2)-C(3) T 3 
C(3) T 3 



hio x 



T 3 \ T 
X ~T 3 ) 7 I r _ 



(22) 



The second line of Eq. (|22|) is a new contribution involving 
the macroscopic value of (No) below the critical temper- 
ature. Since ksT,, ^> hiu x here, it is the dominant contri- 
bution to the fluctuations of the particle number differ- 
ence. It clearly leads to the occurrence of a maximum of 
these fluctuations, at a temperature which remarkably is 
independent of the trap anisotropy: 



-2/3 



0.63. 



(23) 



The corresponding variance is strongly super-poissonian 
in the large atom- number limit: 



[Var(iV L - 7V fl )] r 



(AO 



(24) 

For the parameters of the upper curve in Fig. f2J these 
approximate formulas lead to a maximal variance over 
(A) equal to ~ 27.5, whereas the exact result is ~ 22.2, 
located at T/T c ~ 0.61. We thus see that finite size cor- 
rections remain important even for the large atom num- 
ber (N) = 13000. Fortunately, it is straightforward to 
calculate the next order correction. For the condensate 
atom number we simply expand the summand in (fT4|) up 
to order rj^ , and we recover the known result [l3| : 



(No) 
(A) 



1 - =r 
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T 2 3C(2) fiWrn 

T c 2 2C(3) k B T c 



(25) 



with the arithmetic mean uj m = J^a^a/S. For Coo we 
use the Euler-Mac Laurin summation formula, applied 
to the previously defined function f(x) over the interval 
(1, +00), and we obtain 



21n2 



1 



'h 



(26) 



A^ + o{ n r) 

-2 5 / 2 C(-1/2)/tt ~ 0.374 14]. Using these 
more accurate formulas for (Nq) and Coo in the second 



with A 
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term of the right-hand side of (|12l) leads to an excellent 
agreement with the exact result, see the dashed lines in 
Fig. [2] practically indistinguishable from the solid lines. 
Note that the effect of the —1 correction in (|2"B1) is to 
change the shot noise term 1 in the square brackets (|22|) 
and ((231 into 1 - (N )/(N). 



IV. A PHYSICAL ANALYSIS SINGLING OUT 
THE CONDENSATE MODE 

To investigate the contribution of different physical ef- 
fects on our observable, it is convenient to go back to 
the expression ((3]) and split the field operator into the 
condensate and the non-condensed part: 



ip(r) = (f>{r)a + 6ij>(r) , 



(27) 



where </>(r) is the ground mode wavefunction of the 
harmonic potential. The pair correlation function g( 2 > 
is then expressed as the sum of three contributions, 
g,(2) ( rj r ') = -\-gW -f gfjj, sorted by increasing powers 
of Sip: 



(2) 

9 1 



r <-> r 



2 (r)^ 2 (r')(44aoa o > 
gf} = [0(r)0(r')(aSa o ^ + (r')#(r)) 

+ U 2 (T)(ala 5tp\T')5t()(T')) + ror' 
(<^ + (r)<^V)^( r ')<¥( r )) • 



(2) 

9m 



(28) 

(29) 
(30) 



Averages involving different numbers of operators clq and 
vanish since the system is in a statistical mixture of 
Fock states in the harmonic oscillator eigenbasis. 

(2) 

The term g) ' originates from the condensate mode 
only. Its contribution to V&t(Nl—Nfi) is zero for symme- 
try reasons. This is a crucial advantage, because it makes 
our observable immune to the non-physical fluctuations 
of the number of condensate particles in the grand canon- 
ical ensemble, and legitimates the use of that ensemble. 

(2) 

For the same symmetry reasons, the second line of g IT 

( 2) 

has a zero contribution to Var(A r ^ — Nr). The term g)^ 
originates from the non-condensed gas only. Below T c 
this gas is saturated (z ~ 1) to a number of particles scal- 
ing as T 3 , see (|2"0"|) . and its contribution to V&t(Nl — Nr), 
maximal at T = T c , makes the fluctuations in the par- 
ticle number difference only weakly super-poissonian, as 
already discussed. 

(2) 

Below T c , the first line of g n is thus the important 
term. It originates from a beating between the conden- 
sate and the non-condensed fields. Its contribution to 
Var(A r L — Nr) can be evaluated from Wick's theorem 



Vaxn(N L -N R ) = 4(JV ) 



r)0(r')g (1) (r,r') 



r£i Jr'GL 

iW) 5 (1) (r,r') 



r£L Jr'ER 



(31) 



Using Eq. ([8} and setting z ~ 1 in after some alge- 
bra, we obtain for T < T c 



Vav n (N L - N R ) ~ (N )c o 



(32) 



We can thus give a physical meaning to the mathematical 
splitting (TT2"j) for T <T C : The second term and the sum 
over s in the right-hand side of (|I2[) respectively corre- 
spond to the condensate-non-condensed beating contri- 
bution Va.iii(Ni J — Nr) and to the purely non-condensed 
contribution V&tiu(Nl — Nr). 



V. CLASSICAL FIELD APPROXIMATION AND 
INTERACTING CASE 

In this section we show for the ideal gas that the clas- 
sical field approximation [13, EH, H3] exactly gives 
the high temperature limit (fcgT ^> fvjj x ) of the ampli- 
tude Coo in the condensate-non-condensed beating term 
of Yai{N L - Nr) fl32}. We then use the classical field 
approximation to extend our analysis to the interacting 
case. 



A. Ideal gas: test of the classical field 
approximation 

It is useful to rewrite (|3"Tj) as an integral over the whole 
space introducing the sign function s{x). One then rec- 
ognizes two closure relations on r and r' and obtains: 



\s{x) 



1 



1 



s{x)\(j>) 



(33) 



Correspondingly, in the classical field limit: 

k B T 



^class 



2{<j>\s{x) 



s(x)\<f>). (34) 



Inserting a closure relation on the eigenstates of the har- 
monic oscillator \n) in (|3"3")l . we are then led to calculate 
the matrix elements in one dimension 



x (0\s(x)\n} 1 



2h 

mui x n 



1/2 



^S(0)^-i(0). (35) 



To obtain this result, we introduced the raising opera- 
tor a x of the harmonic oscillator along x and we eval- 
uated the matrix element x (0|[s(a;), a],]\n — l) x in two 
different ways. First, it is equal to n 1 ^ 2 x (0\s(x)\n) x since 
a x \n—l) x — n 1 / 2 !?!}^. Second it can be deduced from the 
commutator [Y(a;),a^] = [H/(2mcj x )] 1 ' 2 S(x) where Y(x) 
is the Heaviside distribution. From the known values of 
0fj(O) (see e.g. Q), we thus obtain: 



3 (2m+l)» 7a! 



mGN 



(2m)! 



2 2m (2m + l)(m!) 



H2 



(36) 
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The equivalent for the classical field is readily computed 
and we obtain: 



nlass 



2 In 2 

Vx 



(37) 



showing that the classical field approximation gives the 
right answer for the dominant contribution to our ob- 
servable. Moreover, going to the first order beyond the 
classical field approximation, that is including the —1/2 
term in the expansion l/[exp(u) — 1] = u^ 1 — 1/2 + 0(u), 
equation (|33[) readily gives the term —1 in (|26[) . In what 
follows we will use the classical field approximation to 
treat the interacting case. 



B. Classical field simulations for the interacting gas 

In Fig|3]we show results of a classical field simulation in 
presence of interactions for two different atom numbers 
(blue circles and black triangles). The non interacting 
case for one atom number (red circles and red curve) is 
shown for comparison. We note that the assumption of 
an ideal Bose gas is nowadays realistic: Recently, the use 
of a Feshbach resonance has allowed to reach a scattering 
length of a = 0.06 Bohr radii [18|. We estimate from 
the Gross-Pitaevskii equation for a pure condensate that 
interactions are negligible if \Ng J |(/)| 4 <C fiw m in where 
w m i„ is the smallest of the three oscillation frequencies 
w a . For the parameters of Fig. [3] (red curve) this results 
in the well-satisfied condition N <C 10 5 . 

In presence of repulsive interactions, the peak of fluc- 
tuations in the particle number difference at T < T c is 
still present, approximately in the same position, but its 
amplitude is strongly suppressed with respect to the ideal 
gas case. Another notable effect is that the dependence 
of the curve on the atom number is almost suppressed in 
the interacting case. 

To compute the normally ordered contribution in 
(j3]) we generate 800 stochastic fields in the canon- 
ical ensemble sampling the Glauber-P function that 
we approximate by the classical distribution P oc 
8 (N - J \ip\ 2 ) exp{-pE[ip,i>*]} where E[tp, ip*] is the 
Gross-Pitaevskii energy functional 



fJ 



m 4 



(38) 



where the coupling constant g = Airfi?a/ra is propor- 
tional to the s-wave scattering length a. The approxi- 
mate Glauber-P function is sampled by a brownian mo- 
tion simulation in imaginary time [2l|. The simula- 
tion results are plotted as a function of T/T£ lms where 
the transition temperature in the classical field simu- 
lations, extracted by diagonalization of the one-body 
density matrix, slightly differs from T c given by (1161) : 
T ciass = L15T 
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FIG. 3: 

for N = 



(Color online) Symbols: Classical field simulations 
17000 (circles) and TV = 6000 (triangles) 87 Rb atoms 
with interactions (blue and black), and without interactions 
for comparison (red). Red solid line: analytical prediction ([9]) 
in the non interacting case. Blue and black solid lines: analyt- 
ical prediction (I54[) in the interacting case. The oscillation fre- 
quencies along x,y,z are tj a /27r[Hz] = (234, 1120, 1473). For 
the interacting case the s-wave scattering length is a — 100.4 
Bohr radii. 



C. Analytical treatment for the interacting gas 

In the interacting case, we perform the same splitting 
as in (|2"T|) except that now the condensate field ipo = 
(iVo) 1 / 2 ^ solves the Gross-Pitaevskii equation 



(hi + gipo)il>o = M) • 



(39) 



For the dominant contribution of to the signal, for 
T < T c , we then obtain 

g$(r,i>) = 2(/V o )0(r)0(r') [(At(r')A(r)) + (A(r')A(r))] 



(JVo) [Ar)(At(r')A(r')} 



r «-> r 



(40) 



where we have neglected fluctuations of Nq, set ao — 
(A ) 1/2 e ie and introduced A(r) = e^Stf;^) [22]. The 
second line in (|40| brings no contribution to the signal 
for symmetry reasons. Note that the terms in </ 2 ' that 
are cubic in the condensate field vanish since (A(r)) = 0. 

In the Bogoliubov [23| and classical field approxima- 
tion, the non-condensed field A has an equilibrium prob- 
ability distribution 

P(A,A*)cxexp|-^ |(A*,A)r?£(^)} (41) 

where the matrix rjC is given in (25j . Splitting A into 
its real and imaginary parts Ar and A/, which turn out 
to be independent random variables, we then obtain the 
probability distribution for A#: 



P(A R ) cx exp -0 / ArHAr 



with 



H = hi + 3gi/jQ — m ■ 



(42) 



(43) 
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The relevant part of g^) , first line in (|40|) , can then be ex- 
pressed in terms of only. Passing through the eigen- 
states of H and proceeding as we have done to obtain 
([3"5 ]) -(|54" ]) . we then obtain 



= 2k B T(<f>\8(x)±a(x)\ 



(44) 



that is the equivalent of (f3~4"f for the interacting gas. We 
then have to solve the equation 



(45) 



Hx = s(x)4> ■ 



In the absence of the sign function the solution 4> a of ((45 
is known 22, 261: 



Mr) = 



d(jVo)TMr) 
(N ) 1/2 n>((No)) 



(46) 



where // is the derivative of the chemical potential with 
respect to (iVo). This can be obtained by taking the 
derivative of the Gross-Pitaevskii equation ([55)1 with re- 
spect to (No). 

In the presence of s(x), the spatially homogeneous case 
may be solved exactly. One finds that the solution x 
differs from s(x)(j) a only in a layer around the plane x = 
of width given by the healing length £. In the trapped 
case, in the Thomas-Fermi limit where the radius of the 
condensate is much larger than £, we reach the same 
conclusion and we use the separation of length scales to 
calculate \ approximately. Setting 



X(r) = /(r)0 Q (r) 



(47) 



we obtain the still exact equation 
h 2 grad <f) a 



2m d> 



grad/ + / = s(a;). (48) 



We expect that / varies rapidly, that is over a length scale 
£, in the direction x only, so that we take A/ ~ d 2 f and 
we neglect the term involving grad /. Also / deviates 
significantly from s(x) only at a distance < £ from x — 0, 
so that <t> a /<t> may be evaluated in x — only. This leads 
to the approximate equation 



Neglecting the deviation of x from s(x)(j) a in the thin 
layer around x = 0, we finally set 



X(r) ~ s(x)</> Q (r) 



and obtain the simple result 



k R T 



(52) 



(53) 



For harmonic trapping where [i scales as {Nq) 2 ^ 5 , and 
assuming that {N ) is well approximated by the ideal 
gas formula (|20| . we then obtain: 



Var (N L - N R ) ~ (N) 



knT r 



m((aO) 



3/5 



T 



(54) 

The analytic prediction (f5~4"]) is plotted as a full line (black 
and blue for two different atom numbers) in Fig[3j We 
note a good agreement with the numerical simulation. 
From (|54p we can extract the position of the maximum 
in the normalized fluctuations: 



T 

as well as their amplitude: 
[Var {N L - N R )} 



980 1 / 3 
14 



0.709 



(N) 



1 + 1.36 



M«iV» 



(55) 



(56) 



Note the very weak dependence of kfjT r / fj,((N)) on the 
atom number, scaling as (N}^ 1 / 15 (27|. We point out 
that our analysis in the interacting case is quite general 
and can be applied to atoms in any even trapping po- 
tential provided that the potential does not introduce a 
length scale smaller than the healing length £. 



D. First quantum corrections to the classical field 

Expanding the non-condensed fields A(r) and A^(r) 
over Bogoliubov modes: 



( A(r) 
^At(r) 



%( r ) 



+ bl 



>\ vj(r) «!(r) 



(57) 



with 



K 2 (y,z) 



d 2 J + f = s(x), 



h k (y,z) _ <f>(0,y,z) 



2m 



§a (0,2/, z) 



20* -U) 



(49) 



(50) 



where U is the trapping potential and the approximation 
of k was obtained in the Thomas-Fermi approximation 
gtpQ ~ fj, — U. Equation ((49|) may be integrated to give 



X(r) ~ s(x)(/)a(r) 



,—k(.V,z)\x\ 



(51) 



and using l|40p . the coefficient Coo giving the dominant 
contribution (f3"2"|) to the signal for T < T c , can be split 



as Cqo 



„th 



With the thermal contribution 



^ = 2^n J |(0| S (a:)(K) + h))| 2 >O (58) 

3 

where fij = (b^bj) = l/[exp(/3ej) — 1], Cj being the en- 
ergy of the Bogoliubov mode j, and the zero temperature 
contribution 

cl = 2j2Ms(x)\v 3 )\ 2 + (4>\8(x)\«i)(vMx)\<l>) ■ (59) 



7 



Performing the classical field approximation that is set- 
ting fij = knT/cj in (|55|) , exactly gives the expression 
(|44| of c^f ss [3l[ . Introducing the first quantum correc- 
tion — 1/2 to the occupation number fij and including the 
quantum contribution, we obtain the first quantum cor- 
rection to c^f ss in the interacting case <5coo = c^, + Sc^ : 

SCoo = ^ ((f>\ S (x) \Vj ) (Vj | S (x) | (j>) - ((f>\ S (x) | Uj } (Uj | s (x) | (j>) . 

(60) 

Using the closure relation JY \ u j)( u j \ ~ l^'X^jl = 1 ~ 
10X01 HI, we obtain 



(61) 



5c n 



This exactly corresponds to the first correction in (|26[) for 
the ideal gas. We conclude that also in the interacting 
case, the first quantum correction to c^f ss has the effect 
of changing the shot noise term (the first term equal to 
one) in the square brackets of (l54l) into I — (N )/(N). 



VI. CONCLUSION 

We have studied fluctuations in the difference of num- 
ber of particles in the left (x < 0) and right (a; > 0) halves 
of a three dimensional harmonically trapped Bose gas as 
a function of the temperature across the critical temper- 
ature for condensation. Both for the ideal gas and the 
interacting gas, fluctuations are weakly super poissonian 
for T > T c . If one lowers the temperature from T c down 
to zero, fluctuations increase, reach a maximum and then 
decrease again as the non-condensed fraction vanishes. 
We have solved this problem analytically for the ideal 
gas case, and we have found an approximate solution in 
the interacting case in the Thomas-Fermi limit when the 
temperature is larger than the quantum of oscillation in 
the trap. Remarkably, the local density approximation 
fails for this problem for the ideal gas. On the contrary, 
we show that the classical field approximation correctly 
gives the high temperature contribution to the fluctua- 
tions in the particle number difference both for the ideal 
gas and the interacting gas. 

For a large ideal gas, the maximum of normalized fluc- 
tuations Var(iV L - N R )/{N) is located at T/T c = 2~ 2 / 3 
independently on the trap oscillation frequencies and its 



amplitude approximately scales as N 1 ^Q/lj x . For the in- 
teracting case in the Thomas- Fermi regime the maximum 
of fluctuations in the relative number of particle subsists 
approximately for the same value of T/T c but its ampli- 
tude is strongly reduced as well as its dependence on N 
scaling as iV" 1 / 15 . 

Finally, we give a physical interpretation of the 
"bump" in fluctuations for T < T c as due to a beating be- 
tween the condensate mode 4>o ( r ) an d the non-condensed 
modes Sip(r), as it is apparent in Eqs. (|29I40I) . For sym- 
metry reasons, only the antisymmetric component of the 
non-condensed field 6iPa(x) = —5iPa(—x) contributes to 
Var(iVi — Nfj). The measurement can then be seen in 

signal 

4> (x) + S\|f (x) 




FIG. 4: The measurement can be seen as a balanced homo- 
dyne detection of the non-condensed field where the conden- 
sate acts as a local oscillator 



a pictorial way as a balanced homodyne detection of the 
non-condensed field where the condensate field acts as a 
local oscillator (Fig|3}. Since one or the other of these 
two fields vanishes for T tending to or T c , the beating 
effect, and thus Var(7Vi — Nr), are obviously maximal at 
some intermediate temperature. 
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